LOCAL SYMMETRY OF HARMONIC SPACES AS DETERMINED BY THE 
SPECTRA OF SMALL GEODESIC SPHERES 



TERESA ARIAS-MARCO AND DOROTHEE SCHUETH 

Abstract. We show that in any harmonic space, the eigenvalue spectra of the Laplace operator 
on small geodesic spheres around a given point determine the norm |Vii| of the covariant derivative 
of the Riemannian curvature tensor in that point. In particular, the spectra of small geodesic 
spheres in a harmonic space determine whether the space is locally symmetric. For the proof 
we use the first few heat invariants and consider certain coefficients in the radial power series 
expansions of the curvature invariants \R\^ and |Ricp of the geodesic spheres. Moreover, we 
obtain analogous results for geodesic balls with either Dirichlet or Neumann boundary conditions. 
We also comment on the relevance of these results to constructions of Z.I. Szabo. 



1. Introduction 

For a compact closed Riemannian manifold S the spectrum of S is the eigenvalue spectrum, 
including multiplicities, of the associated (positive semi-definite) Laplace operator A acting on 
smooth functions. A central question of inverse spectral geometry asks to which extent the 
geometry of S is determined by its spectrum. The so-called heat invariants ak{S) of S are 
examples of geometric invariants which are determined by the spectrum of S; indeed, they are the 
coefficients in the famous asymptotic expansion by Minakshisundaram-Pleijel, 

oo 

TV(exp(-tA)) ~ (47rt)-^*'"(^)/2^afc(5)t^ 

A:=0 

for t 4, 0. The first few of these coefficients are given by 



ao(S) = vol(S), aUS) = - [ scaldvols, aafS) = — / (Sscal^ - 2|RicP + 2[i?|- 

6 Js 360 Js 



')dvols , 



where seal, Ric, and R denote the scalar curvature, the Ricci operator, and the Riemannian 
curvature operator of S, respectively. In general, each ak{S) is the integral over S of certain 
curvature invariants; see [S] for more information. 

Nevertheless, there exist many examples of pairs or families of isospectral Riemannian manifolds 
(i.e., sharing the same spectrum) which are not isometric, sometimes not even locally isometric; 
see, for example, the survey article [lOj. Still, many questions remain open; for example, it is not 
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known whether a locally symmetric compact closed Riemannian manifold can be isospectral to a 
locally nonsymmetric Riemannian manifold. 

On the other hand, the geometry of geodesic spheres plays an interesting role in Riemannian 
geometry. Chen and Vanhecke [4] formulated the following general question: To what extent 
do the properties of small geodesic spheres determine the Riemannian geometry of the ambient 
space? For example, Gray and Vanhecke [12] studied the information contained in the volume 
function of small geodesic spheres and investigated the question whether a Riemannian manifold 
whose geodesic spheres have the same volumes as spheres in euclidean space must necessarily be 
flat (answering this question in the positive under various choices of additional assumptions). 

In the context of inverse spectral geometry, an interesting special version of the above question 
is: To what extent do the spectra of small geodesic spheres in a (possibly noncompact) Riemannian 
manifold M determine the geometry of M ? For example. Theorem 6.18 in [4J uses the information 
contained in the heat invariants qq and ai of small geodesic spheres (viewed as functions of the 
radius) and concludes local isometry of manifolds with adapted holonomy to certain model spaces 
under the assumption that all small geodesic spheres around each point are isospectral to the 
corresponding geodesic spheres in those model spaces. 

In order to arrive at such and similar results, one uses radial power series expansions of curvature 
invariants, both of the ambient space and of the geodesic spheres. In general, even the first few 
coefficients of such expansions become very complicated; see, for example, the various formulas 
in [12] or [1]. One setting in which a quite restrictive geometric assumption on the ambient space 
makes the calculations considerably easier is the setting of harmonic ambient spaces. 

A manifold is called harmonic if the volume density function of the geodesic exponential map 
is radial around each point. The notion of harmonicity was first introduced by Copson and 
Ruse [5] and intensively studied by Lichnerowicz [13]; see also [16]. Chapter 6 of the book by 
Besse f2] gives a useful survey of properties of harmonic spaces. One of the important facts 
about harmonic spaces is that they are Einstein [2j and hence analytic [7J (the latter result was 
not yet known when [2j was written). A locally symmetric manifold is harmonic if and only if 
it is flat or of rank one; the famous Lichnerowicz conjecture postulated that, conversely, each 
harmonic space is locally symmetric; i.e., satisfies Vi? = (this condition is classically known to 
be equivalent to the condition that the local geodesic symmetries around each point be isometries). 
For the case of compact manifolds with finite fundamental group the Lichnerowicz conjecture was 
proved by Szabo [TS]; however, Damek and Ricci gave examples of noncompact homogeneous 
harmonic manifolds which are not locally symmetric in infinitely many dimensions greater or 
equal to seven [6]. These spaces are usually referred to as Damek-Ricci spaces; see for more 
information. 

Specializing the above question about the information contained in the spectra of small geodesic 
spheres to the setting of harmonic spaces, we are able to prove in the present paper that the spectra 
of small geodesic spheres in a harmonic space determine whether the space is locally symmetric 
(see Corollarv 11.21 below). More precisely, we obtain: 

Main Theorem 1.1. Let Mi and M2 be harmonic spaces, and let pi G Mi, p2 € M2. If there 
exists e > such that for each r € (0,e) the geodesic spheres Sr{pi) and Sr{p2) are isospectral, 
then \ VR\l^ = \VR\l^. 

Corollary 1.2. Let Mi and M2 be harmonic spaces. Assume that the hypothesis of Theorem \l.l\ 
is satisfied for each pair of points pi € Mi, p2 G M2. Then Mi is locally symmetric if and only if 
M2 is locally symmetric. 
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In particular, note that in the case of locally homogeneous harmonic spaces Mi and M2, the 
hypothesis of Theorem 11.11 imphes that Mi is locahy symmetric if and only if M2 is locally 
symmetric. Actually, all known examples of harmonic spaces are locally homogeneous; it is an 
open question whether there exist harmonic spaces which are not locally homogeneous. 

Interestingly, our result implies that certain pairs of geodesic spheres which were claimed to be 
isospectral by Szabo in [20], [21] are actually not isospectral. In fact, Szabo considered (as the 
featured examples in a more general construction) geodesic spheres in certain symmetric spaces Mi 
(namely, quaternionic hyperbolic space of real dimension 4m > 12) and in certain associated 
locally nonsymmetric Damek-Ricci spaces M2 of the same dimension (see also Remark 12. 4p . He 
stated that every pair of geodesic spheres Sr{pi) C Mi and Sr{p2) C M2 of the same radius was 
isospectral. Since these ambient manifolds Mi and M2 are harmonic and homogeneous, and Mi 
is locally symmetric while M2 is not. Corollary 11.21 immediately implies that Szabo's result cannot 
be correct. Note that it was Fiirstenau |8J who first discovered that actually there was a gap in 
Szabo's isospectrality argument. The question of whether that proof could be repaired or not had 
since remained open; our result settles this question in the negative. 

The incorrect examples of geodesic spheres mentioned above had the notable property that one 
is homogeneous and the other not. While it remains unknown whether a homogeneous metric on a 
sphere can be isospectral to a non-homogeneous one, Szabo in an earlier article [T9] did construct 
a pair of isospectral metrics, only one of which is homogeneous, on the product of a sphere and a 
torus (those results are not affected by the error in the later papers). 

We obtain analogs of our above results for geodesic balls endowed with either Dirichlet or 
Neumann boundary conditions; see Theorem 15.11 and Corollary 15. 2[ Similarly as above, this 
implies that Szabo's examples in [20], [21] of isospectral geodesic balls (of any given radius) 
in quaternionic hyperbolic space of real dimension at least 12 and in certain associated locally 
nonsymmetric Damek-Ricci spaces were erroneous. 

Note that nevertheless there do exist isospectral pairs and even continuous families of isospectral 
metrics on spheres and balls; the first such examples were due to Gordon [11) . 

In order to prove Theorem 11.11 we use the heat invariants ao{Sr{p)) and 02 (5,. (p)) of geodesic 
spheres in harmonic spaces. In particular, we study the coefficients of in the radial power series 
expansions of ^^i^g^^ fs^^p) \Ric^ l"^ dvols^^p^ and ^^i^g;^ /^^(p) IR^]"^ dvols^^p), where Ric^ and 
denote the Ricci operator and the Riemannian curvature operator of Sr{p)- From the form of these 
coefficients (see Proposition 13. 21 and its mean value version Proposition 14. 2p . we are able to conlude 
that the heat invariants oq and 02 of Srip), viewed as functions of r, together determine the value 
of |Vi?p at the midpoint p. Note that the same is not true for aQ{Sr{p)) alone; see Remark 12.41 
Moreover, in the harmonic setting, the function r ai{Sr{p)) = 5/5 (^p) scal^ dvols^(p) does 
actually not contain more information than r 1-^ ao{Sr{p)) = vol{Sr{p)); see Remark 12.31 So it is 
indeed necessary for our purpose to consider a2{Sr{p))- 

Our computations rely heavily on the harmonicity of the ambient space. Note that they are 
related to certain more general computations in [12] and [1]; for example, Theorem 8.1 of [1] 
actually includes a kind of analog to our Proposition 13.21 and this even for general, not only for 
harmonic manifolds; however, that theorem contains information only on the coefficients of 
with j < 0, while we need the coefficients of r"^. In fact, in the harmonic case, the lower order 
coefficients turn out to be determined already by the function r i-> aQ{Sr{p)) = vol{Sr{p)); see 
Proposition 13.21 and Remark 12. 3( i). 

This paper is organized as follows: 
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In Section [2] we gather the necessary background on harmonic spaces, mostly following [2j; 
in particular, we recall Ledger's recursion formula for the power series expansion of the second 
fundamental form of geodesic spheres, and the resulting curvature identities in harmonic spaces. 

In Section [31 we study the coefficients of for j < 2 in iRic'^''^^'' Ig-^^p^^^^ and l-^'^''^'^'' lexp(rM) 
for unit tangent vectors u of harmonic spaces, using the power series expansion of the second 
fundamental form and its radial covariant derivative, as well as the Taylor series expansion of the 
Riemannian curvature tensor. Proposition 13.21 is the main result of this section. 

Section H] is devoted to the proof of the Main Theorem 11.11 In preparation for this, we first 
derive a mean value version of Proposition 13. 2| see Proposition 14.21 

Finally, in Section \5\ we prove the analog of Theorem 11.11 for geodesic balls. We consider the 
heat coefficients of geodesic balls in harmonic spaces and show that the functions r i— )• ao{Br{p)) 
and r i— )• a2{Br{p)) (either for Dirichlet or for Neumann boundary conditions) together determine 
the value of |Vi?p at the midpoint p of the balls. More precisely, we show that the coefficient of 
in the radial power series expansion of the quotient a2{Br{p))/aQ{Br{p)) is a sum of a nonzero 
multiple of |Vi?|p and of terms determined by the function r i— )• aQ{Bj.{p)). 



2.1. Volume density and the shape operator of geodesic spheres. 

In the following, let M be a complete, connected, n-dimensional Riemannian manifold. For p G M, 
let expp = exp \ ^ : TpM — > M denote the associated geodesic exponential map. For a vector 
V G TpM we denote by 7„ the geodesic with initial velocity v. Identifying Ty{TpM) with TpM, 
we regard the differential (i(expp)^ as a linear map from TpM to T^xp^M. We denote parallel 
translation along 7^ by P^^* : T^^(^g^M r^^(j)M. Given any unit vector u G 5'i(0p) := {u G 
TpM I \u\ = 1} and r G M, we consider the volume density 



with w J-U. Thus, for each r G {0,i{p)), where i{p) denotes the injectivity radius of M at p, 



is the infinitesimal volume distortion at u of the map 

Si{Op) 3 u^ju{r) = exp{ru) G Sr{p), 

where Sr{p) C M denotes the geodesic sphere of radius r around p. Let (7u(r) denote the shape 
operator of Sr{p) at exp(rn); that is. 



where V is the Levi-Civita connection of M and denotes the outward pointing unit normal 



2. Preliminaries 




(1) 



Vu{r) := r^-^Ouir) 
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where R is the Riemannian curvature tensor of M, given by R{x, y)z = —'Vx'^yZ+'Vy'VxZ+'V\^x,y]^: 
and where R^ := .)^- (Note that here we use the same sign for R as Besse f^.) Let 

Cu{r) := rau{r). This endomorphism field along 7M|(oj(p)) smoothly extendable to r = by 
Cu(0) := lu-, where is defined by Iu{u) = and /u||^|X = Idj^jx. Moreover, from ([3]) one can 
derive Ledger's recursion formula for the covariant derivatives of at r = (see, e.g., |4j): 

{k - i)c(t\o) = -Hk - - ^ (^)cW(o)cr^)(o) 

(k) 

for all A; G N, where R^ is the k-th covariant derivative of the endomorphism field R-y.^ along 7^ 

at r = 0. This formula allows one to successively compute the ci!'\o) in terms of the endomor- 

(k) 

phisms Ru of TpM. Forming the Taylor series of Cu and dividing by r, one obtains (see, e.g., 

p;f o au{r) o p^^; = - \Ru - jR'u - i^K + ^RuRuy 

1 „,„ 1 „ „, 1 



(4) " ^36^^" + 72^"^" + ^KRuy 

— f Ro, H RiiR'J, H R!!,Rii -\- R[,R[. -\- Ro,RiiRo,^T^ 

468 " 210 " 210 " 112 " '' 945 ^ 

+ 0^). 

2.2. Curvature identities in harmonic spaces. 

The manifold M is called a harmonic space if for every p £ M the above function 6u does not 
depend on it G 5i(0p). An equivalent condition is that for all r € (0,z(p)), the geodesic spheres 
Sr{p) have constant mean curvature (recall equations ([l]), ([2])). For more information on harmonic 
spaces see |16j or [2]. If M is harmonic then the function 6^ does in fact not even depend on p\ 
that is, there exists 9 : [0, 00) — M such that 

9u{r) = e{r) 

for all u G TM with \u\ = 1. Moreover, even the local or infinitesimal versions of the above 
condition imply that the manifold is Einstein p] and therefore analytic [7j. Hence, the local or 
infinitesimal versions of the above conditions are equivalent to the global versions. Since 6u{r) 
depends only on r, so does Vu{r) and hence Tr((TM(r)). From this one can successively derive, 
using the expansion (j4]): 

Proposition 2.1 (see [2], Chapter 6). If M is harmonic then there exist constants C,H,L G M 
such that for all p ^ M and all u G TpM with \u\ = 1: 

(i) Tr(i?u) = C ; in particular: 

(fi) Tr(i?l''^) = id for all k en. 

(iii) Ti{RuRu) = H; in particular: 

(iv) TtIRuR'J = and 

(v) TriRM = -Tr(KK). 

(vi) Tr{32RuRuRu - 9R'M = L. 

In fact, taking traces in Q, one has in the harmonic case: 

(5) Tr(f7„(r)) = in - 1)- - -Cr - —Hr^ 3— Lr^ + O(r^) 
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for r 4 and all u G TM with \u\ = 1. Note that Proposition I2.1( i) just says that the Einstein 
constant of M is C; that is, Ric = CId on each TpM. Recall that the Ricci operator is defined 
by (Ric{x),y) = Tr(R(x, .)y) for all x,y TpM and all p S M. From Proposition 12.11 one can 
further derive: 

Proposition 2.2 (see [2], Chapter 6). If M is harmonic, then for the above constants C,H,L 
and each p €z M: 

(i) {R{x, . ) . , R{y, ■) ■) = |((^ + 2)-fr — C'^){x, y) for all x,y ^ TpM ; in particular. 

(ii) \R\l = lni{n + 2)H-C'). 

(iii) 32(nC3 + ^C\R\l + lR{p) - R{p)) - 27\VR\l = n{n + 2)(n + 4)L. 

^ o 

Here, the functions R,R(z C°°{M) are certain curvature invariants of order six which are 
defined as follows: If {ei, . . . , e„} is an orthonormal basis of TpM and Rijki '■= {R{ei, ej)ek, ee), 
then 

^ip) •= ^ RijkiRjaebRaibk, R{p) '■= ^ RijkiRkeabRabij ■ 
i,j,k,£,a,b i,j,k,£,a,b 

Note that the term nC^ in Proposition [22f iii) reads nC"^ in the corresponding equation 6.67 in [2], 
but this was obviously a misprint (note that curvature terms of different order cannot occur here) ; 
see also formula (3.1) in [22]. 

Proposition 12. 2( iii) will be used in Section U together with the following formula which actually 
holds in any Einstein manifold; see formula (6-7) in jl4] or formula (11.3) in [12]: 

(6) -^A{\Rf) = 2C\Rf - R-AR+IVR]"^, 

where A denotes the Laplace operator on functions, that is. A/ = — '^i{ei{eif) — {'Veiei)f) for 
local orthonormal frames {ei, . . . , Cn}- (Again, there is a misprint in two of the coefficients in the 
corresponding formula 6.65 in [2j.) If M is harmonic, then the left hand side of (0) is zero by 
Proposition I2.2( ii) . Finally, we recall the following well-known observations which will be used in 
Section H) 

Remark 2.3. Let M be an n-dimensional harmonic space with volume density function 9 as 
above. 

(i) For any p S M, the volume of the geodesic sphere Sr{p) with < r < i{p) equals the 
volume ujn-i of the standard unit sphere S^~^ in M" multiplied by the factor 

v{r) := r"-i^(r) 

Note that v{r) = Vu{r) for each unit vector u € TM, where Vu is the function defined in ([T]). 
The function v determines the volume growth function v' /v of the geodesic spheres, and thus it 
determines, by ([2]), the function Tr(iTu(r)) (which is independent of u). By ([5]), the function which 
associates to small values of r the volume of geodesic spheres of radius r in a given harmonic 
space M determines the constants C, L (and of course n) associated with M. 

(ii) Let seal = nC denote the scalar curvature of M. Let p G M, fix some r G (0, i(p)), and 
let scal'^ denote the scalar curvature function of Sr{p)- A routine calculation using the Gauss 
equation shows that for each unit vector u € TpM we have 

scal^(exp(™)) = seal - 2(Ric(7«(r)), 7«(r)) + (Tr(a„(r)))2 - Tr(a„(r)2) 
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which by the Einstein condition and equations ^ and ([3]) impUes 

scal^(exp(r'u)) = (n - 2)C + {v'{r)/v{r)f + Tr«(r)) + Tr(i?^^(^)) 

= (n - 2)C + {v'{r)/v{r)f + (v'/vYir) + C = {n - l)C + v"{r)/v{r). 

Therefore, geodesic spheres in M have constant scalar curvature, and the respective constant 
depends only on the radius, not on the midpoint. Finally, using (i) one concludes that the 
function which associates to small values of r the scalar curvature of geodesic spheres of radius r 
is determined already by the function which associates to small values of r the volume of geodesic 
spheres of radius r. 

Remark 2.4. As mentioned in the Introduction, the aim of this paper is to show that in harmonic 
spaces, the heat invariants ao{Sr{p)) = vol{Sr{p)) and a2{Sr{p)), viewed as functions of r, together 
determine | Vi?|p . This is not the case for oq alone, as manifested by certain pairs of Damek-Ricci 
spaces. A Damek-Ricci space AN is a certain type of solvable Lie groups with left invariant 
metric, namely, the standard 1-dimensional solvable extension of a simply connected Riemannian 
nilmanifold N of Heisenberg type. The volume density function of is radial and depends 
only on the dimensions of N and its center [6] ; see also the book [1] . Within the class of Damek- 
Ricci spaces, there exist pairs of symmetric spaces AN and locally nonsymmetric spaces AN' 
where A^ and N' have the same dimension and so do their centers. (In fact, certain such pairs 
AN and ^A^' were the ambient manifolds used by Szabo in |20] . [21j : recall the Introduction.) 
In particular, geodesic spheres of the same radius in AN and AN' have the same volume. This 
shows that in harmonic spaces, the function r i— ?> vol{Sr{p)) alone does not determine |V-R|p . In 
turn, Remar k |2 . 3 1 shows that in harmonic spaces, the function r i— )• ao(5'r(p)) = vol{Sr{p)) already 
determines the function r i-^ ai{Sr{p)) = ^ Js (p) scal'^ dvolg^^j,-^ . Therefore, we need to consider 
0'2{Sr{p))- The next section gives some necessary preparations for this. 

3. Radial expansions of |Ricp and for geodesic spheres in harmonic spaces 

In this section we will describe a certain coefficient in the radial power series expansions of the 
curvature invariants |Ricp and of geodesic spheres in harmonic spaces. First we need the 
following lemma. 

Lemma 3.1. Let M he an n- dimensional harmonic space, and let C and H he the constants 
from Proposition \2.1\ Let p £ M, and let S := Sr{p) he a geodesic sphere around p with radius 
r G (0,z(|?)), endowed with the induced Riemannian metric. Let u he a unit vector in TpM , let 
a := o"t((r) he as in Section\^ and write a' := cr^(?'). Let and Ric"^ denote the curvature 
tensor, resp. the Ricci operator, of S. Then in the point q := exp(ru) ^ S we have: 

(i) iRiC^I^ = (n - l)C^ + 2C7(TV(a))2 + {Tv{a)fTT{a^) + 2CTr(a') + 2TY(a)Tr((JCj') + Tr(cjV'), 

2 " 

(ii) = (n - 4) ((n + 2)H - C^) + 4H + 2{Tv{a^)f - 2TY(a^) + 4 ^ TV(a o i?(ei, . )aei) , 

i=i 

where {ei, . . . , e^} is an orthonormal basis of TqM. 

Proof, (i) Let v be the outward pointing radial unit vector field as in Section [2l From the Gauss 
equation one easily derives the following formula whose analog is valid for submanifolds of codi- 
mension one in arbitrary Riemannian manifolds: 

Ric^ = (Ric - Ry^ + Tr(cj)(T - o-^)|t,s 
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Using the Einstein condition and the Riccati equation ([3]), this formula becomes in our situation: 

Ric^ = (CId + Tr(a)a + cT')|T,5 

(see also [H], p. 67). Now one obtains the desired formula immediately, keeping in mind that 
both a and a' are symmetric and annihilate Vq . 

(ii) Choose an orthonormal basis {ei,...,e„} of TqM such that ei = Vq. For all i,j,k,i £ 
{2, . . . , n} we have by the Gauss equation (recall our sign convention for R): 

{R^{ei,ej)ek,ee) = {R{ei,ej)ek,ee) + {aei,ek){crej,ee) - {aej,ek){crei,ei). 

Squaring both sides and forming the sum over i,j,k,£, while recalling that a is symmetric and 
annihilates ei, we get 

n 

l^^lg= Yl (^(ei,ej)efc,e^)^ + klVP + kPkP 
i,j,k,e=2 

n n 

- 2\a^f + 2 X] (^(ei,ej)crej,crej) - 2 ^ {R{ei,ej)aej,aei) 

i,j=l i,j=l 
n n 

= Yl («(ei,ej)efc,e^)2 + 2(Tr(c72))2_2Tr((j^) + 4^Tr((joi?(ei, OcTCi). 

i,j,k/=2 i=l 

The desired formula now follows from the fact that the first sum on the right hand side is equal 
to — 4|ii(ei, .).p+4|i?eiP which bv Proposition 12. 2l fi). (ii) and Proposition 12. iT iii) becomes 
l{n-4){{n + 2)H -C^) +AH. □ 

Using the radial power series expansion of a together with the previous lemma, we will make 
conclusions concerning the first few coefficients of the radial expansions of |Ric^|2 and {R^l"^. The 
following proposition will be the key of the proof of the Main Theorem 11.11 Actually we will use 
only the statements about a2 and in this proposition. 

Proposition 3.2. Let M be an n- dimensional harmonic space, and let C, H, and L be the 

constants from Proposition \2.1[ Let p E M , and let u be a unit vector in TpM. Then 

iRic^'^^^lcxpCru) = a-4r"^ + a_2r"^ + oq + a2iu)r'^ + O(r^) and 

iR'^-^^^llMru) = /5-4^~' + /3-2r-2 + f3o + f32{uy + 0(r3) 

for r J, 0, where the coefficients ai and (3i for i € {—4, —2, 0} are constants depending only on n, 
C , and H. Moreover, 

a2(u) = 0(2 + ■^Ti(R'R') and 
lb 

4 " 



/32(w) = /32 + -XlTr(i?«oi?(ei, .)Rue^), 



i=l 



where a2 and (32 are constants depending only on n, C, H , and L, and where {ei, . . . , e„} is an 
orthonormal basis of TpM . 

Proof. We use Lemma [3T] together with the power series expansions ([5]) of a := cJu(r) and 
Trfcr). Let us first consider |Ric'^''*-^^|^ and the individual contributions of the nonconstant 
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terms in Lemma l3.1( i) to its expansion. By ([5]) we have 

(Tr(a))2 = ((n - l)i - \cr - ^Fr^ - -^Lr^f + 0{r^) 

for r J, 0. Moreover, from the expansion @ and Proposition 12.11 one gets 

(7) Tr(CT2) = (n - 1)1 - -C + —Hr"^ + -^Lr^ + Oir^). 

w V ; V V2 3 15 3024 ^ ^ 

Further, 

Tr(a') = -^Tr(c7) = -(n - 1)1 - 1(7 - l^r^ + 0{t') 
ar 6 15 

2Tr(c7a') = ^Tr(cj2) = -2(n - 1)1 + ^Hr + -^Lr^ + O(r^) 
ar r'^ 15 756 

by ([7]). Using these expansions and ([5]), one easily checks that each of the the four individual 
terms 2C(Tr(o"))^, (Tr((T))'^Tr(cT^), 2CTr(o"'), and 2Tr((T)Tr((Tcr') appearing on the right hand side 
of Lemma [3.1( i) has the property that the corresponding coefficients of r~'^, r"^, depend only on 
n, C, the coefficient of depends only on n, C, H, L, and the coefficients of r^^,r~^,r vanish. 
It remains to consider the term Ti{a'a') in Lemma l3.1( i). From (jl]) we get 

o a' o =-!/.- ^Ru - iK - + Yf^RuRuV 



by dSD, and 



( R^f^ H RiiR[', H R[',Rii H R[,Ri, H R^iRjiR^i^r^ 

^168 ' 42 " « ' 42 « " ' 112 " " 189 ^ 



+ 0(r5) 

and thereby, using Proposition I2.lt 



Tiia'a') = {n-l)^ + ~ + ^H 

3r2 45 



The coefficient of r in the latter expansion is 



Tt(R'R') + Tt(RuRuRu) 

840 ^ " 945 ^ " 



which by Proposition I2.1l fvi) turns out to be 



62 ,37 9 • 62 , , , , , 31 1 , , , , 

-L + ( + )Tt(R'R') = L + —Tr(R'R'). 



32-945 ^840 32-945^ ' " 15120 16 

This concludes the proof of the statements concerning the expansion of |R'ic'^''^^^|gxp(-j,„) • 

We now turn to and study the individual contributions of the nonconstant terms 

in Lemma I3.1( ii) to its expansion. Squaring ([7j), we see that in the expansion of the term 
2(Tr(cj2))2 the coefficients of /--^ j,-2^^o depend only on n,C,H, the coefficient of depends 
only on n, C, H, L, and the coefficients of r~^,r~^,r vanish. 
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Regarding the term — 2Tr((T^) we obtain from 

p-;° ° o p^: = - - + {-Ik + 1^5) + {-Ik + 1(r.k + kkj)-- 

for r 4, 0. Using Proposition 12. II we get 



-2Tr(a^) = - 2(n - 1)1 + - 

3 45 



The coefficient of in the latter expansion is 

10 1 28 

420 ^ " 945 ^ " " 

which by Proposition I2.1( vi) equals 

It remains to consider the term 4^^-^Tr(cr o i?(e.j, ■)(Tei) in Lemma l3.1l fii). We make some 
preliminary observations. For k € Nq, let R^^\ resp. Ric^'^^ denote the k-ih. covariant derivative 
of the curvature tensor, resp. the Ricci operator, along 7^ at r = 0. We will use the the Taylor 
series expansion of the Riemannian curvature tensor along 7„ (recall that M is analytic): 

00 ^ 

(9) p;^>R,.ir)op',: = Y.'yR^'^ 



k=0 



Moreover, Ric^'^^ = for k > 1 since M is Einstein. Note that Ric| j^j^^ = X^iLi similarly 
on each T^^^j.^M if we extend {ei,... ,en} parallelly along 7^. For any € No we have, using 
Proposition 12 II 



5^ TY(/„ o i?«(ei, . )Iuei) = Tr(/„ o Ric('=)) - Tr(/„ o i?('=)) 



(10) 



i=l 



Moreover, 



(11) 



Tr(Ric('=)) - (Ric('^)n, u) - Tr(i?W) = ^ " J; 



^T¥(i?„oi?W(e,, .)/„e,) =Tr(i2,oRicW) -TV(i?„i?(f)) 



1=1 



C^-H, k = 0, 

-TriRuRi^^), k>l, 
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^Tr«oi?W(e,, .)4e,) =Tr(i?:,oRicW)-TrK4^)) 



(12) 



0, k = 0, 

-Tt{B!^b!u^), k>l, 



(13) ° • )Iuei) = Tt{RI o Ric) - Tt{RIR^) = + Tr(i?;0, 

4 = 1 

n 

(14) ^^{RuRu o ii(e„ . )/„ei) = Tr(i?„i?„ o Ric) - Ti{RuRuRu) = CH- Tr(i?„i?„i?„). 

i=l 

Note that for any pair of symmetric endomorphisms F, G of TpM we have 

n n 

(15) j;Tr(Foi2W(e„ OGe,) =^Tr(Goi?W(e„ . )i^e,) 

i=l 1=1 

by the symmetries of the curvature operator. Keeping the expansions and ([9]) in mind, we see 
that the expression in (jlOp contributes only to the coefficient of in the expansion of Y2i=i Tr(co 
R{ei, . )aei) , the expression in ([TT]) contributes to the coefficients of r° and (and higher order), 
the expressions in (jl2p . ()13p . ()14p contribute to the coefficient of (and higher order). The 
only additional contribution to the coefficient of is given by the sum of Tr(i?u o i2(ej, . )i?„ej). 
Recalling (|15p (and multiplying i?'^'^) by 1/A;!), we obtain from @, I2.ir iv). ([9]), and the above 
observations: 



4f;Tr(^7oi?(e„ . )ae,) = 4((n - 2)-^ - ^(C^ - /7) 
1=1 



+ 



LS"^^'^^"^") + " ^)Tr(KK) - ^CH + ATr(i?,ii„i?„) 
1 " 



i=l 



By Proposition I2.1( v). the coefficient of r in the latter expansion is 

-^CH - ATr(«) + ^Tr{R^RA) + ^ V Tr(i?„ o R{e,, . )i?„e,) 
45 15 45 9 ^-^ 

i=l 

By Proposition I2.1( vi). the two terms involving Tv{R'^R'^) and Tv{RuRuRu) become 
" •^+(-ll + ^)Tr(«) = T^^-^Tr«i?:,). 



32-45 ' 15 32-45' ' " 180 12 

Combining this with the result for the r^-coefficient of —2Tr (cr'*) from ([8]), we conclude that the 
terms involving Tr(i?^i?^) in the coefficient of in the power series expansion of |-^?'^''^^^lexp(r«) 
cancel each other, and the only remaining term apart from those which depend solely on n, C, H, L 
is I Ya=i Tr(i??i o R{ei, . )Ruei), as claimed. □ 
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Remark 3.3. For the purpose of the proof of the Main Theorem II .11 in Section [U which we will 
perform using the heat invariants ao{Sr{p)) = vol{Sr{p)) and a2{Sr{p)) = ggg js (p)(^(^*^^^'^)^ ~ 
2|Ric'^P + 2[i?'^P)duoZ5^,(p), we would actually not have needed the exact statement of the previous 
proposition - which might, however, be interesting in its own right. Rather, we could have 
restricted our attention to the term Th:{a'a') in the expression of jRic^j^ ^.^^ in Lemma b-lf i). and 

to the last two terms in the expression of ^^^^ in Lemma l3.1( ii). In fact, even without the 

explicit calculation of the expansion of the other terms, one easily sees that those are determined 
by the volume function r ^ ao{Sr{p)) = vol{Sr{p)) of the geodesic spheres (which is just the 
function v multiplied by the volume of the standard unit sphere, see Remark 14.11 below). More 
precisely, in the spirit of Remark 12.31 we obtain 

2C{Tr{a)f = 2C{v'/vf, 
{Tr{a)fTr{a') = {v'/vfi-{v'/vy-C), 
2CTi{a') = 2CTr(a)' = 2C{v'/vy, 
2Tr(a)Tr(aa') = 2v' /v ■ ^(Tr(cj2))' = v' /v ■ {-{v'/v)"), 
2{Ti{a^)f = 2{-{v' /v)' -Cf. 



4. Proof of the Main Theorem 

In this section we will first derive an integrated version of Proposition 13.21 Using this and the 
heat invariants ao, 01,02 of geodesic spheres in harmonic spaces we will then prove our Main 
Theorem 1 We need the following general remark on mean values. 

Remark 4.1. In any harmonic space M, the average (or mean value) of a smooth function / 
on a geodesic sphere Sr{p) (with < r < i{p)) is the same as the average of /(exp(r.)) over 
the unit sphere 51(0^) in TpM. More explicitly: Let con-i denote the volume of the (n — 1)- 
dimensional standard sphere. In particular, w„_i is the volume of Si{Op). Recall from Section [2] 
that 9{r) = 6u{t) is independent of u (and even of p) by harmonicity. We have 

\o\{Sr{p)) = r'^'^e{r)ujri-i = v{r)ujn-i , 
and for any smooth function / on Sr{p), 



f dvols^(p) = —— / f{exp{ru))v{r)du 

/(exp(rn)) du. 



V0l{Sr{p)) Jsrip) v{r)uJn--l J Si{Op) 

1 



^n-l JSi(Op) 

Now we can give an "integrated" version of Proposition 

Proposition 4.2. Let M be an n-dimensional harmonic space, and let C, H, and L he the 
constants from Proposition {2711 Let p ^ M . Then 

— , / \Ric^''^'''^\'^dvoL(r,) = a-4r~^ + Q!_2r~^ + ao + 02?^^ + 0(r^) and 

V0l{Srip)) Js,{p) 

/ Ii?^'-(^')|2rf^oZ5.(p) = /3-4r-^ + /3-2r-2 + /3o+^2^' + 0(r3) 
vol(5r(p)) Js^(p) 
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for r 10, where the coefficients ai and f3i for i G {—4, —2, 0} are the constants from Proposition \3.S\ 
depending only on n, C , and H . Moreover, 

- 3 , ,2 

ao = Q!2 H 77 r Vi? L and 

16n(n + 2)(n + 4)' 

''^ = ft + MS^I^'^l- 

where 6i2 and (^2 o,re constants depending only on n, C , H , and L. 

Proof. For any unit vector u in TpM, let a2{u) and I32{u) be the coefficients from Proposition 13.21 
Using that proposition and Remark 14.11 we only need to show that 

_ I f I f 1 

7x2 '■= / a2{u) du = a.2-\ / — Tr(i?^i?^) du and 

^2 ■■= [ /?2(u) du = P2 + [ ^ V Tr(i?„ o R{ei, . )R^ei) du 

are of the claimed form, where 02, /32 are as in Proposition 13.21 For 02 this follows immediately 
(with 0(2 ■= 0(2) from the following formula (see the proof of Theorem 5.7 of [15j; details of the 
computation can be found on p. 170 of [12j): 



(16) / MKK 

JSiiOp) 



du = — z — " -iVRl'i 

n(n + 2)(n + 4)' 



This confirms the statement concerning 02 • 

We now consider (32- Writing u = J27=i "^i^i ^"^^ ^ijke = {R{ei, ej)eki ^i) we have 



n 



'^Tr{RuO R{ei, .)Ruei) = ^ {R{ei,ej)ek,ei){Ruei,ek){R 

i=l i,j,kl=l 

(17) 



n n 

1^ ^ RijklRaibkRcjdl UaUbUcUd- 
a,b,c,d=l i,j,k/=l 



Note that the integral oi UaUhUcUd over 5i(0p) is zero whenever {a,b,c,d} contains at least three 
different elements. Abbreviating Aabcd ■= J2ijke=i^ijkiRaibkRcjdi we have, using the Einstein 

condition and recalling the definition of R and R from Section [2) 

n n n n 

^aabb = ^ RijkiRaiakRbjbl = ^ Rijki^ikSjl = ^ Rijij = nC^ , 

a,b=l a,b,i,j,k,£=l i,j,k,£=l *ij=l 

n n n 

A-abab = ^ RijkiRaibkRajbl = RijklRaibkRjalb = R{'P)i 

a,fe=l a,b,i,j,k,£=l a,b,i,j,k,£=l 

n n n ^ 

y^ ^abba = ^ RijklRaibkRbjal = RijklRaibkRjbia = R{p) — -^R-ip), 

a, 6=1 a,b,i,j,k,£=l a,b,i,j,k,£=l 

where for the last equality we have used formula (2.7)(vi) of |17j : see also formula (2.15) of [12]. 
Let C i?" denote the (n — l)-dimensional standard sphere. Note that Jg„^i u^u^ du = ^^^'^j 
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and fgn-1 uf du = '^^1^2) • From the above equations and (fT7|l we thus obtain 



/I ft 

/ V'Tr(_R„ o R{ei, . )Ruei) du 

= / , i^aabb + -^abab + -^abba 

a,be{i,...,n} "^-f a=l 

n 

E [^-^^ + ^abab + Aabba]^^^^ = [uC^ + 2R{p) - -R{p)) 



■ 

a,b=l 

Hence, 

(18) a = A + + 2fl(p) - ifl(p));^ = ft + ^ + - ifl(p));^ 

Recall from Proposition I2.2r ii) . (iii) and equation ([6|) that 

ll2R{p) — 32R{p) = 27|Vi?[p + some constant depending only on n, C, H, L and 
R{p) + 4:R{p) = |Vi2|p + some constant depending only on n,C,H, 
using which one easily computes that 

— -R{p) + -R{p) = -|Vi?|p + some constant depending only on n,C,H,L. 
9 9 8 

Thus we conclude from ([18]): 

where /32 is a constant depending only on n,C, H, L. □ 
Proof of the Main Theorem ll.lt 

Let Ml, M2 be harmonic spaces, pi G Mi, p2 G M2, and assume there exists e in the interval 
(0, min{i(pi), i(p2)}) such that for each < r < e the geodesic spheres Sr{pi) and Sr{p2) are 
isospectral. Then dim Mi =: n = dimM2, and the heat invariants of the geodesic spheres coincide: 



akiSriPi)) = ak{Srip2)) 

Ipi I * -"-1^2 ■ 



for each r G (0,e) and all k G Nq. We want to deduce that |Vi?|p = \\7R\3, . Actually this will 



follow using just oq and 02- 

Reformulating the problem, let M be an n-dimensional harmonic space and p M. We want 
to show that for any e G {0,i{p)), the two functions 

^Pk : (0,e) Br^ak{Sr{p)) GM 

with k G {0, 2} together determine the value of |Vi?|p . By Remark I2.3f i). the function 

ipo -.r ^ ao{Sr{p)) = vol{Sr{p)) = v{r)ujn-i 

determines the constants C, H, L associated with M (see Section [2]). Recall that the scalar curva- 
ture scal"^ =: scal'^'' of Sr{p) is constant on the manifold Sr{p), and that the function u : (0, e) ^ M 
determines, by Remark I2.3( ii). the function (0, e) 9 r 1— )• scal'^'' G K. In particular, the function 
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(/?o = v'j-'n-i also determines the function (0,e) 3 r ^^^^{scal^ dvols^(p) = ipoir)-{scal'^'^)'^ 

By 



Mr) = a2{Sr{p)) = (5(scal^)2 - 2|Ric^p + 2\R^\^) dvols^^^) 

-aOU JSr{p) 

it follows that the functions <y9o and 992 together determine the function 

(0,e) 9 r ^ I I - |Ric^|2)d^o/5,(p) G M. 

By Proposition 14.21 the r^-coefficient in the power series expansion of this function is the sum of 
the term 

1 3 2n + 5 i„^|2 



8n(n + 2) 16n(n + 2)(n + 4) / ' 16n(n + 2)(n + 4) 

and /32 — Q!2- Recall that the latter is a constant depending only on n,C, H, L, and is thus 
determined by ipo. We conclude that the functions ipo and ip2 together determine |V-R|p , as 
claimed. □ 



5. Geodesic balls 

In this section we will prove the following version of the Main Theorem 11.11 for geodesic balls: 

Theorem 5.1. Let Mi and M2 he harmonic spaces, and let pi S Mi, p2 € M2. If there exists 
e > such that for each r G (0,e) the geodesic halls Br{pi) and Br{p2) are Dirichlet isospectral, 
then \S/R\l^ = iVi?!^^. The same holds if the assumption of Dirichlet isospectrality is replaced by 
the assumption of Neumann isospectrality. 

This theorem implies the corresponding analog of our Main Corollary 11.21 

Corollary 5.2. Let Mi and M2 he harmonic spaces. Assume that the Dirichlet isospectrality 
hypothesis of Theorem I5.il is satisfied for each pair of points pi G Mi, p2 S M2. Then Mi is 
locally symmetric if and only if M2 is locally symmetric. The same holds if the assumption of 
Dirichlet isospectrality is replaced hy the assumption of Neumann isospectrality. 

For the proof of Theorem 15. II we will use the heat invariants for manifolds with boundary. Let 
M be an n-dimensional Riemannian manifold, and let i? C M be a compact domain with smooth 
boundary. If A denotes the Laplace operator on B with Dirichlet boundary conditions then there 
is an asymptotic expansion 

Tr(exp(-tA)) ~ (47rt)-«/2 af 

k 

for t 4, 0, where k = 0,0.5, 1, 1.5, . . . ranges over the nonnegative half integers (see p]). For the 
Laplace operator on B with Neumann boundary conditions the analog of this formula holds with 
certain coefficients (B). The coefficients (B) (resp. {B)) are given by certain curvature in- 
tegrals over B and dB. One has a^{B) = {B) = vol(5) and 0^5(5) = -0^5(5) = -^vol{dB) 
(see [3]). In the proof of Theorem 15.11 we will use the explicit formulas for a^(i?) and (B) 
from [3]. Let v denote the outward pointing unit vector field on the boundary dB of B, and let 
a = Viy be the associated shape operator. Let seal, Ric, R always refer to the usual objects on M 
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(not to the ones associated with the induced metric on dB). Then 
360 



/ (-12A(scal) + 5scal^ - 2|Ric|^ + 2|i?p) dvols 
Jb 



+ / (l8zv(scal) + 20scal • Tr(o-) - 4Tr(i?^)Tr(o-) + 12Tt{R^ o a) 
JdB ^ 



40 88 320 

4Tr((Ric - R,) o a) + -{TT{a)f - -Tr{a)Tv{a') + -^Tv{a') ) dvoks 



a^(B) = — I (-12A(scal) + 5scal2 -2|RicP + 2|i?HduoZB 
360 [Jb 

+ / (-42zy(scal) + 20scal • Tr(o-) - 4Tr(i?^)Tr(o-) + 12Ti-(fi^ o a) 

- 4Tr((Ric - R^) o a) + ^{Tr{a)f + 8TY(CT)Tr(CT2) + Y^r{a^)) dvoks 

If M is harmonic then, by the results of Section [21 the previous formulas simplify to 

«?(^) = 1^ Uol{B) ■ {5inCf - 2nC^ + ^n((n + 2)H - C^)) 
360 I 3 

(19) + / (20nCTr{a) - 8CTr{a) + 16Tr(i?^ o a) 

+ ^(Tr(-))^ - f Tr(<T)Tr(.2) + ^-^Tr{a')) dvoles 

a^{B) = ^ [vol(i?) ■ (5(nC)2 - 2nC^ + ^n((n + 2)H - C^)) 

(20) + / f 20nCTr(f7) - 8CTr(CT) + 16Tr(i?^ o a) 

+ y (Tr(a))3 + 8Tr(a)Tr(a2) + f Tr(a3)) dvolsB 

In the proof of Theorem 15.11 we will follow a similar strategy as in the proof of Theorem 11.11 To 
this end, we need some preliminary results in the special case that B = Brip) with r € (0,«(p)) 
and M is harmonic. We remark - without going into details this time - that one can compute in 
this case, using equation Proposition 12. H and the Taylor series expansion 



OO I, 



k=0 



that the r^-coefficient in the power series expansion of Tr(i?,^ocr) equals —j^L-\--^Ti{R'^R'^), and 
that the r^-coefficient in the power series expansion of Tr(o"^) equals 35^^ — ^Tr(i?^i?^). (That 
the contributions of Tr(i?^i?^) in these terms are negatives of each other can also be checked as 
follows: Using ([3]) twice, we have Tr(i?^oCT) + Tr(o-3) = -Tr(o-V) = -^Tv{a'^y = ^TT{Ru + a'y = 
^Tv{a)" whose r'^-coefficient indeed depends only on L by ([5]).) Using p^ . we conclude that the 
r^-coefficient in the power series expansion of r H- ^q^^g (p^) fs (p) Tr(i?i, o a)dvolg^(p-^ is 



(21) 



1 



1 



1440 32n(n + 2)(n + 4) 



\VR\l, 
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Similarly, the r^-coefficient in the power series expansion of r i— )• ^q/^g (p)) fsr{p) '^^i^^)'^^^hr{p) 

^^^^ 30240^ " 32n(n + 2)(n + 4) '^^'^ ' 

Proof of Theorem I5.lt 

Let Ml, M2 be harmonic spaces, pi S Mi, p2 G M2, and assume there exists e in the interval 
(0, min{i(pi), i(p2)}) such that for each < r < e the geodesic spheres Br{pi) and -Brfe) are 
Dirichlet isospectral (resp. Neumann isospectral). Then dim Mi =: n = dimM2, and the heat 
invariants of the geodesic spheres coincide: 

af (i?,(pi)) = af (S,(p2)), resp. {Br{pi)) = (i?.^)) 

for each r G (0,e) and all k G Nq. We want to deduce that |Vi?|p^ = |Vi?|p2- Actually this will 
follow using just oq and 02- (We remark without proof here that, viewed as functions of r, the heat 
coefficients ao.5, ai, and ai.5 do actually not contain more information than ao in our situation.) 

We first consider the case of Dirichlet conditions. Similarly as in the proof of Theorem 1 1.1^ we 
reformulate the problem as follows: Let M be an n-dimensional harmonic space and p G M. We 
want to show that for any e G (0,i(p)), the two functions 

: (0,6) 9r^af(i?,(p)) GM 
with k G {0, 2} together determine the value of |Vi?|p . Note that the function 

V'? : r ^ a^{Br{p)) = vol(S,(p)) 
determines its own derivative which is just 

r ^ \o\{Sr{p)) = v{r)LOn-l 

(see the previous section). By Remark l2.3( i). we conclude that ifj^ again determines the constants 
C, L associated with M. Moreover, the function : (0, e) — )• M determines the radial functions 
Tr(o-) = v'/v and Tr:{a'^) = -{v' /v)' - C (compare Remark [33]). By dUD it now follows that V'o' 
and V'2' together determine the function 

1 I' 320 

Recalling (|2ip and (j22p . we see that the r^-coefficient in the power series expansion of the latter 
function is the sum of 

1 , 320. 2 _ 1 |^p|2 

32n(n + 2)(n + 4)^^^ ~ ^j|V^lp - 42n(n + 2)(n + 4)' 

and a term depending only on L. Since L is determined by tpQ , we conclude that tpQ and V'2' 
together determine |Vi?|p , as claimed. 

In the Neumann case, letting (r) := {Br{p)), we again have V'o^ (r) = Yol{Br{p)) = ipQ{r). 
Proceeding exactly as in the Dirichlet case, using (f20]l this time, we see that ipQ and ip2 together 
determine the sum of 



32n(n + 2)(n + 4)^ 3^' 6n(n + 2)(n + 4) ' 

and a term depending only on L. Hence they determine |Vi2|p , as claimed. □ 
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